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ISOMETRIC IMMERSIONS OF SURFACES WITH TWO CLASSES OF 
METRICS AND NEGATIVE GAUSS CURVATURE 

WENTAO CAO, FEIMIN HUANG, AND DEHUA WANG 


Abstract. The isometric immersion of two-dimensional Riemannian manifolds or sur¬ 
faces with negative Gauss curvature into the three-dimensional Euclidean space is studied 
in this paper. The global weak solutions to the Gauss-Godazzi equations with large data 
in L°° are obtained through the vanishing viscosity method and the compensated com¬ 
pactness framework. The L°° uniform estimate and H~^ compactness are established 
through a transformation of state variables and construction of proper invariant regions 
for two types of given metrics including the catenoid type and the helicoid type. The 
global weak solutions in L°° to the Gauss-Godazzi equations yield the isometric 
immersions of surfaces with the given metrics. 


1. Introduction 

Isometric embedding or immersion into of two-dimensional Riemannian manifolds 
(or surfaces) is a well known classical problem in differential geometry (cf. Cartan 
in 1927, Codazzi [12] in 1860, Janet [30| in 1926, Mainardi [3l| in 1856, Peterson [H] in 
1853), with emerging applications in shell theory, computer graphics, biological leaf growth 
and protein folding in biology and so on (cf. [211139] ). The classical surface theory indicates 
that for the given metric, the isometric embedding or immersion can be realized if the first 
fundamental form and the second fundamental form satisfy the Gauss-Codazzi equations 
(cf. |5l [35l|36l E]) • There have been many results for the isometric embedding of surfaces 
with positive Gauss curvature, which can be studied by solving an elliptic problem of the 
Darboux equation or the Gauss-Godazzi equations; see m and the references therein. 
When the Gauss curvature is negative or changes signs, there are only a few studies in 
literature. The case where the Gauss curvature is negative can be formulated into a 
hyperbolic problem of nonlinear partial differential equations, and the case in which the 
Gauss curvature changes signs becomes solving nonlinear partial differential equations of 
mixed elliptic-hyperbolic type. Han in [26| obtained local isometric embedding of surfaces 
with Gauss curvature changing sign cleanly. Hong in [5H] proved the isometric immersion 
in of completely negative curved surfaces with the negative Gauss curvature decaying 
at a certain rate in the time-like direction, and the norm of initial data is small so 
that he can obtain the smooth solution. Recently, Ghen-Slemrod-Wang in [8| developed 
a general method, which combines a fluid dynamic formulation of conservation laws for 
the Gauss-Codazzi system with a compensated compactness framework, to realize the 
isometric immersions in with negative Gauss curvature. Christoforou in HI] obtained 
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the small BV solution to the Gauss-Codazzi system with the same catenoid type metric as 
in [8]. See [THl [lU [20l [28l [311 HU SH [50] for other related results on surface embeddings. 
For the higher dimensional isometric embeddings we refer the readers to |2l[3l[l[9l[ini[22l 
[251 [Ml [Ml SOI 03] and the references therein. 

Recall that in Chen-Slemrod-Wang [8|, the isometric immersion of surfaces was 
obtained for the catenoid type metric: 

ds^ = E{y)dx^ + cE{y)dy^ 

and negative Gauss curvature 

K{y) = -koE{y)~^^ 

with index /3 > ^/2, c = 1 and > 0. The goal of this paper is to study the isometric 
immersion of surfaces with negative Gauss curvature and more general metrics. To this 
end, we shall apply directly the artificial vanishing viscosity method (IE]) and introduce 
a transformation of state variables (see e.g. [27]). The vanishing viscosity limit will be 
obtained through the compensated compactness method ([Il[7l[Ill[l5l[l6l[23l[29l[32l[33l 
EZl SI]). We shall establish the estimate of the viscous approximate solutions by 
studying the Riemann invariants to find the invariant region (see e.g m) in the new state 
variables. Then we prove the H~^ compactness of the viscous approximate solutions, and 
finally apply the compensated compactness framework in [8] to obtain the weak solution 
of the Gauss-Codazzi equations, and thus realize the isometric immersion of surfaces into 
M^. In the new state variables, using the vanishing artificial viscosity method, we are able 
to obtain the weak solution of the Gauss-Codazzi system for two classes of metrics, that 
is, for the catenoid type metric with /3 > \/2 and c > 0, and for the helicoid type metric 

ds^ = E{y)dx^ -|- dy^ 

and 

K{y) = -koE{yf 

with a < —2 and /cq > 0- An important step to achieve this development is to find the 
invariant regions of the viscous approximate solutions for the wider classes of metrics, for 
which the introduction of new state variables (tt, v) plays a crucial role. We note that the 
L°° solution of the Gauss-Codazzi equations for the given metric in yields the 
isometric immersion from the fundamental theorem of surfaces by Mardare [351136] . 

The paper is organized as follows. In Section 2, we introduce a new formulation of the 
Gauss-Codazzi system and provide the viscous approximate solutions. In Section 3, we 
establish the estimate for the two types of metrics to get the global existence of the 
equations with viscous terms. In Section 4, we prove the compactness. In Section 
5, combining the above two estimates and the compensated compactness framework we 
state and prove our theorems on the existence of weak solutions for the surfaces with the 
catenoid type and helicoid type metrics. 


2. Reformulation of the Gauss-Codazzi System 

As in m and |8] , the isometric embedding problem for the two-dimensional Riemannian 
manifolds (or surfaces) in can be formulated through the Gauss-Codazzi system using 
the fundamental theorem of surface theory (cf. Mardare [351 [36] ). 


ISOMETRIC IMMERSIONS OF SURFACES 


3 


Let ri C be an open set and {x,y) G Q. For a two-dimensional surface defined on n 
with the given metric in the first fundamental form: 

I = Edx^ -|- 2Fdxdy + Gdy"^, 


where E, E, G are differentiable functions of (x, y) in n. If the second fundamental form 
of the surface is 

E = hiidx^ -|- 2hi2dxdy + h22dy'^, 

where hii,hi 2 = /i 2 i )^22 are also functions of {x,y) in Q, then the Gauss-Codazzi system 
has the following form: 

Ly = rl^L - 2Tl^M + rj^N, 

< N:,-My = -Tl2L + 2T\2M-T\^N, ( 2 . 1 ) 

LN = K. 

\ 

Here and in the rest of the paper Dy stand for the partial derivatives of corresponding 
function □ with respect to x,y respectively. In (12.ip . 


L = 


/ill 


M = 


h 


12 




N = 


h22 

VW\' 


\g\ =EG-F^>0, 


K = K{x,y) is the Gauss curvature which is determined by E,F,G according to the 
Gauss’s Theorem Egregium ([Ulllb]), F^^- {i,j,k = 1,2) are the Christoffel symbols given 
by the following formulas (EZl): 


F 

F 

F 


1 

11 


1 

12 


GEx — 2FFx -\- FEy 

2{EG - F2) 
GE„ - FGx 


— F^ — 

— i 21 — 


2{EG - F2) ’ 


1 

22 


2GFy - GGx - FGx 
2{EG - F2) 



2EFx — FEy — FEx 

2{EG - F2) ’ 



FGx — FEy 
2{EG - F2) ’ 

^ - 2FFy + FG 
2iEG - F2) 


X 


( 2 . 2 ) 


As in Mardare [351136] and Ghen-Slemrod-Wang [ 8 ], the fundamental theorem of surface 
theory holds when (hij) or L,M,N are in for the given E,F,G G C^’^, and the 
immersion surface is Therefore it suffices to find the solutions L,M,N in L°° to the 
Gauss-Codazzi system to realize the surface with the given metric. 

In this paper, we consider the isometric immersion into of a two-dimensional Rie- 
manian manifold with negative Gauss curvature. We write the negative Gauss curvature 
as 

K = - 7 ^ (2.3) 

with 7 > 0 and 7 G C'^’^(fl), and rescale L, M, N as 


~ L ~ M ~ N 
L = -, M = —, N=—. 

Ill 

Then the third equation (or the Gauss equation) of system (12.ip becomes 


(2.4) 


= - 1 . 


( 2 . 5 ) 
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The other two equations (or the Codazzi equations) of system (|2.ip become 

M^-Ly = tl^L - 2fl^M + 


where 


^2 

22 


Nx- 

II 

pi 

1 22 

i 2f 12 

M - 

-fliiV 

1 

— ^ 22 


■p2 

1 12 

= r?2 + 


- 

11 - 

= r?i, 


7 



27 



— J- 22) 

pi 

1 12 

— 1 15 

+ 1 ]L 

pi 

111 

1—11—1 

II 

i^. 




27 

7 


( 2 . 6 ) 


(2.7) 


Consider the following viscous approximation of system (I2.5I) - (I2.6I1 with artificial vis¬ 
cosity: 


Ly — Mx — eLxx — ^22^ “I" 2f ^2-^ ~ I 

My — Nx = eMxx + 1^22-^ “ 2 T\ 2 M + f 
LN -M"^ = -1, 


( 2 . 8 ) 


where e > 0. Our goal is to apply the vanishing viscosity method to the smooth solutions 
of (j2.8p to obtain the L°° solution of ()2.5p - (l2.6p . The eigenvalues of the system ()2.6p for 
L / 0 are 

^ -M ± 1 
A± = 


L 


and the right eigenvectors are 
A direct calculation shows 


r± = (1,-A±)'^ 


-Mil -1 -Mil 

VA± • 'f'± =-^— • 1 H—^ ^— = 0 , 

^ -L2 L -L 


thus we can take \± as the Riemann invariants. 
Introduce the new variables: 

M 


u := —- 


L 




then 


Thus 


L = 


M = 


—u 


N = 


V? - 


T — __i 
— 2 ’ 


~ 2vl - VVx 

■L/T.T. - o 


Nx = 


{2uux — 2vvx)v — {vr — v^)vx 2uvux — v^Vx — u^Vx 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


and 


“ „2 ’ 


Mx = 


UVx — VUx 


Mxx = 


UVxx — VUx 


2uVx — 2vUxVx 


( 2 . 12 ) 
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Substituting (12.111) and (12.121) into the system (|2.8I) . we get 


—Vy uVfr — vu^, 2evi — evv, 


XX f^2 


■pz 

1 oo — 1 O — -L 1 


9 2 

U -oU— V 


and 


22 “ ^-^ 12 “ -^11 

V V V 


UVy — VUy 2uVUx — v'^Vx — U^Vx 

-- 

_ euvxx - £VUXX 2euvl - 2evUxVx 1 , oui ^ ri 

r 22 “ + 12 77 + 11 y ■ 


y2 y3 

Multiplying (I2.13P by —we get 


2ev, 


Vy + UVx — VUx = £Vxx - - + r22'*^ + 2ri2UV + — V‘^)v 


and multiplying ()2.14p by — u, one has 


UVy 


+ Uy + 2uUx — VVx — 


U^Vx 


euv. 


— Sl/'.r.T 


V 


2sUxVx 2euVx -i ^- 1,0 9n 

—^ - r^2 - 2ri2« - _ ^2^ 


Substituting (I2.15p into (12.161) . we obtain 


Uy UUx VVx —£'^xx ^22 ^^12^ ^11 ^ ) 

+ t22U + 2f^2^^ + — v'^)u 

2sVrU 


- 


- fL + (fi 2 - 2 f} 2 )ix + ( 2 f ?2 - fhy 


I T'l 2 I t '2 ( 2 2 \ 

Therefore we have the following system in the variables {u,v): 

2£VxUx 


Uy + {uUx - VVx) =£Uxx “ 


“ ^22 + (f 22 — 2t\2)u 


+ (2f'i2 ~ + fiiU^ + Tii{u'^ — v‘^)u, 

2r7;2 

X I T '2 I ot 2 I t 2 / 2 2 \ 

- ^22^ “ 1 “ ^llv^ — ^ 


Vy + {uVx - VUx) =£Vxx - 

V 

We note that the L°° estimate for the solutions L, M and iV to (|2.8p with L 
equivalent to the L°° estimate for the solutions u and v to (|2.17l) with u > 0. Set 

j(u, V) = -tl2 + (fi2 - 2f 12)11 + {2tl2 - f }i)n2 + fliu2 + f2,(u2 - V^)u, 

g{u, v) = f| 2 U + 2f I 2 W + - v'^)v, 

then the system (|2.17p becomes 


Uy + {uUx - VVx) = fiu, v) + euxx - 


2eUxVx 


Vy + {uVx - VUx) = giu, v) + evxx - 


2£V. 


(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

> 0 is 

(2.18) 

(2.19) 
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The local existence of (12.81) is standard, and the global existence can be proved if the 
L°° boundedness of L, M and N (or equivalently the L°° boundedness of u and v) is 
established. The L°° uniform bound will be established in the next Section [3l The global 
existence of solution (L, M) to (12. 8 p is equivalent to the global existence of solution (u, v) 
to (j2.19l) which will be proved in Section [5l 


3. Uniform Estimate 

In order to establish the L°° bound of u and v, we need to derive the equations of the 
Riemann invariants of (12.191) . First we rewrite the system (12.191) in the following form: 


tiy 


U 

— V 


1 

_ 1 


1 

_ 1 


^ ^XX y 

1 

_ 1 

+ 

— V 

u 


1 

_ 1 


1 

1 _ 

+ 

2ev'^ 

evxx -^ 


The eigenvalues of (13.ip are 

Xi = u — V, X2 = u + V, 

and the Riemann invariants are 


W = U + V, z = u — V. 

Multiply (13.Ij) by {wu,Wv) or {zu,Zy) to obtain the equations satisfied by the Riemann 
invariants: 

2eVxWx j./ \ ! \ 

Wy + XiWx = ewxx - J[u,v) +g{u,v), 

o ^ (3-2) 

zevxZx 

Zy + X2ZX = ezxx - \- j{u,v) - g[u,v). 

V 

Then at the critical points of w, the first equation of (j3.2h becomes 

swxx + f{u, v) + g{u, v) = 0 , 

and at the critical points of 2 :, the second equation of (13.21) becomes 

szxx + f{u, v) - g{u, v) = 0 . 

Hence, by the parabolic maximum principle (see mm). we have 

(a) w has no internal maximum when f{u,v) + g{u,v) < 0 , 

(b) w has no internal minimum when f{u,v) + g{u,v) > 0, 

(c) 2 ; has no internal maximum when /(tt, v) — g{u, v) < 0 , 

(d) 2 ; has no internal minimum when f(u, v) — g{u, v) > 0 . 

In order to find the invariant region of w, z, we need to analyze the source terms in ()3.2p . 
that is, the signs of f{u, v) + g{u, v) and f{u, v) — g{u, v). We shall consider two types of 
surfaces that have special metrics with F = 0. 
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3.1. Catenoid type surfaces: G{y) = cE{y),F = 0, c > 0. For the surfaces with 
metrics of the form: 

G{y) = cE{y), F = 0 (3.3) 

and (EZI), we can easily calculate 


with constant c > 0, from the formulas of F^^- in 
that 


F^ 

22 2E' 

^22 = 0 ) 


^12 — 0 , 

F^ 


i 11 — 


-E' 

rli = 0 , 


and thus 


E' 


7^2 _ _|_ 7 

-*■ 22 — ^ ^ W 


■fs _ n 

12 — U, 


F^ — 
^ 11 — 


-E' 

2 ^' 


^22 ~ 0 ) 


2 E;^ 27 ’ 


fli = 0 . 


If we assume 

where a is constant, then, by (12.3p . 


i E' 

— = 

7 E 


„ E' K' 


that 


IS 


K[y) = -koEiy)^^, 


with some constant ko > 0. From the expressions of f{u,v), g{u,v) in (|2.18|) . one has 


(3.4) 


and then 


, . E' 1 E' 2 2\ 

g{u, v) = —{a + -)v - ^{u - v )v. 


E' 

f{u, v) + g{u, v) = {cu - c{2a + l)i; + (u^ - v‘^){u + v)) , 

E' 

f{u, v) - g{u, v) = {cu + c{2a + l)v + (u^ - v^){u - v)) . 


Set 


(pi(u, v) = cu — c{2a + l)u + (tt^ — v‘^){u + v), 
(P 2 {u, v) = cu + c{2a + l)u + (rt^ — v‘^){u — v). 
In particular, when a = —1, 

(pi(u,v) = (u^ - + c)(u + v), 

(P2(u, v) = {u^ — + c){u — v). 

If we assume 

E' 


(3.5) 


(3.6) 


(3.7) 
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Figure 1. Graphs of level sets of w and 2; 


then the signs of 

f{u,v)+g{u,v), f{u,v)-g{u,v) 

depend only on ipi{u,v),ip 2 {u,v) respectively. 

We now derive and sketch the invariant regions. 

When a = —1, first we can easily sketch the level sets of w,z in the u — v plane in 
Figure 1. To obtain the invariant region we need to sketch the graphs of 

(piiu,v) = 0, (P 2 iu,v) = 0 

in the u — v plane in Figure 2. Now let us find the invariant region in the upper half-plane. 
We draw a straight line parallel to tt — u = 0 passing through the point 

C = {0,5) 

with 0 < h < ^/c intersecting the hyperbola at the point 

A = {uo,uo + 6), 

and similarly, we get the point 

B = {-uo,uo + 6), 

where 

uo = uo(-l) = (3-9) 

We then draw a straight line perpendicular to tt — u = 0 passing through the point A and 
another straight line perpendicular to tt + u = 0 through B. Then the two lines intersect 
the u—axis at the same point 

c — 
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Figure 2. Graphs of (pi{u,v) = ^,ip 2 {u,v) = 0 when a = —1 



Figure 3. Invariant region when a = — 1 


see Figure 3. We see that the square ACBD in Figure 3 is an invariant region. Therefore 
we get the estimate of {u,v), that is, 


c-52 


<u< 




6 <v < 

0 


25 


25 


(3.10) 
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Figure 4. Graphs of ipi{u,v) = Q,^p 2 {u,v) = 0 when a < —1. 


When a < — 1, the graphs of v) = 0, n) = 0 in the u — v plane look like the 
curves marked with 1 and 2 respectively in Figure 4. Similarly to the case a = —1, we 
draw a straight line parallel to u — u = 0 passing through the point 


C = {Q,5) 

with 0 < (5 < \J—2ca — c intersecting the curve (pi {u, u) = 0 at the point 

A = {uo,uo + 6). 

Then we draw a straight line parallel to u + u = 0 passing through the point C = (0,5) 
intersecting the curve ip 2 {u,v) = 0 at the point 

B = {-uo,uo + 5). 

Moreover, we draw a straight line perpendicular to u — u = 0 passing through point A and 
another straight line perpendicular to u + u = 0 passing through point B. Then the two 
lines intersect the u—axis at the same point 


D — (0,2tto + 5). 


Thus, the square ACBD in Figure [5] is an invariant region, and now the L°° estimate of 
(u, v) is 

— uq < u < uq, 6 < V < 2uo + 5, (3T1) 


where 

0 < 5 < V—2ca — c, 
ca + 25^ — V (Ao? — 4ca5^ — 4c5^ 


Uq = uo{a) 


45 


(3.12) 
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Figure 5. Invariant region when a < — 1 


Example 3.1. For the following special catenoid type surfaces with 

1 


F(y) = (ccosh 0))^ , G{y) = 


c2(/32-1)2 


E{y), 


and 


Kiy) = -c\l3^ - l)E{y)-f^ , 

where c / 0, /3 > \/2 are constants, one has F'(y)' > 0 whenever y > 0, and -E(y)' < 0 
whenever y < 0. All the conditions (13.4p and ()3.8I) for the above invariant region are sat¬ 
isfied when y < 0. If we take = {(x, y) : x € M, —yo < y < 0}, where yo > 0 is arbitrary, 
then in the equations (|2.19p are parabolic for y—time like. Therefore, when the initial 
value (u(x, —yo), x(x, —yo)) is in the square ACBD, the parabolic maximum/minimum 
principle ensures that the square ACBD is an invariant region, which yields the L°° esti¬ 
mate. We notice that the surface is just the classical catenoid when /3 = \/2- Indeed, for 
the special catenoid type metric, the surface is given by the following function: 

r(x,y) = (ri(x,y),r 2 (x,y),r 3 (x,y)). 


with 


ri(x,y) = (^ccosh 
r2(x,y) = (^ccosh 


1 


sin(x), 

cos(x). 


r^ix.y) = 


f 


(^2 _ 1)2 


ccosh I - 


4 - 2 / 3 '^ 

B'^-1 


dt. 
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3.2. Helicoid type surfaces: E{y) = B{y)‘^,F = 0,G{y) = 1. For the surfaces with the 
metric of the form: 


E{y) = B{y)\ E = 0, G(y) = 1, B{y) > 0, 
we can also calculate that 


2 9 9 -E' 1 1 E' 

r 22 = 0 ) ^12 = 0 ) ~ ” 2 ^’ ^22 = 0 ) ^12 = ^ 


and 


— f) 

^ 22 ~ —’ ^ 12 ~ '-'j 


7 


^ 11 — 


-E' 


^22 ~ 0 ) 


2E ^ 27 


f 11 = 0 . 


Then 


f/ ) ^ ( 2 2\ 

f{u,v) = -—u-—{u -v)u, 


7 


E' 


Assuming 


g{u, v) = —V - —{u^ - v^)v. 
7 ^ 


y B' 

~ ~ “"b ’ 


7 


and using E{y) = B{y)‘^ and (|2.3I) . we have 


that is, 


with some constant ko > 0. Thus, 


E' K' 


K{y) = -koE{y)^, 


2B' 


and 


T 1 =0. 


) 11 


Now set 


f{u, v) = -^u - BB'{u^ - y )n, 

g{u, v) = a—v — BB'{u^ — v‘^)v, 

B 


‘2B^ B^ 

f{u, v) + g{u, v) = - —u + a—v - BB'{u^ - v‘^){u + v), 

B B 

2tB^ B^ 

f{u, v) - g{u, v) = ——u - a—v - BB'{u^ - v‘^){u - v). 


u = Bu, V = Bv, 

w = Bu + Bv = u + V, z = Bu — Bv = u — v. 


(3.13) 


(3.14) 


(3.15) 
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Wy + XiWx = ewxx - 


2evxWx 


Note that 


~\~ X 2 ZX — 


B' 


V 

2eVxZx 

V 


+ (/(^^, v) + g{u, v))B + —w 


B 


B' 


(3.16) 


+ {f{u, v) - g{u, v))B + —z. 


f{u, v) + g{u, v) = {2Bu - aBv + {{Bu)'^ - {Bv)‘^){Bu + Bv)) , 


then 


B^ 

Ri{u, v) :={f{u, v) + g{u, v))B + —w 

B^ B' 

= - — {2u - av + {u^ - v‘^){u + h)) + — (n + v) 

= - ^ {u - {a + l)v + {u^ - v‘^){u + v)) , 


and similarly, 


B^ 

R 2 {u, v) :={f{u, v) - g(u, v))B + —z 
B^ 

= - — (it + (a + l)h + {u^ - v^){u - h)) . 

Since B[y) >0 is given, it remains to find the invariant region of w, z. As in the catenoid 
case in Subsection EH we need to analyze the signs of iii and i? 2 - If we assume 

B' < 0, 

or equivalently, from E = B^, 

E' < 0, (3.17) 

and set 

v) := u — {a + l)v + {u^ — v^)(u + v), 

V’ 2 (it)h) := it + (a + \)v + {u^ — v^){u — v), 

then the signs of Ri{u,v), R 2 {u,v) depend on the signs of -01 (it, h), 1 /^ 2 (it, h) respectively. 
Now we end up with a situation similar to the catenoid case in Subsection 13.11 with c = 
1,0 = 2a, and the corresponding ipi (i = 1, 2) in Subsection 13.11 are 

Lpi{u, V) = ijjiiu, V), ip2{u, V) = ljj2{u, V). 

We can find the invariant regions just as in catenoid case in Subsection 13.11 Indeed, the 
invariant region of (it, h) looks like the same as the invariant region ACBD in Figure El 
when o = —2 if we replace the u — v plane by the u — v plane and uq by ito = uo{—2) 
defined below; and looks like the same as the invariant region ACBD in Figure El when 
a < —2 in the u — v plane, and thus we omit the sketch of the invariant regions. From the 
invariant region, when a = —2, 

1 


1-52 1 _ 

- < U < -;-. 

- - 2 (f ■ 


i' 


2(1 


(3.18) 
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and when a < —2, 


— uq < u < uq, S < V < 2uo + 5, 


where 

0 < (5 < V—a — 1, 

—a — 46^ + — 8ad^ — 165^ 

uo = uo{a) = - — -. 

Note that 0 < B(y) G we can easily obtain the L°° boundedness of u,v. 


(3.19) 


(3.20) 


Example 3.2. For the helicoid surface with 

E(y) = c^ + y^, F(y) = 0, G(y) = 1, K(y) = - f 

[c^ + y^)^ 

where c 7 ^ 0 , we see that 

B(y) = \/c^ + y^, 

and B(yy > 0 for y > 0, and B(yy < 0 for y < 0. If we take 

^ = {(a:, y) : X eR, -yo <y <0}, 

where yo > 0 is an arbitrary constant, then in Q, the equations in (I2.19jl are parabolic 
for y—time like. Therefore, when the initial value (n(a:, —yo), u(x, —yo)) is in the square 
ACBD, we have the invariant region for the solutions. We note that the function for the 
helicoid in is r(x,y) = (y sinx, y cos x, ex). 

Remark 3.1. We can see from Figures [3] and [5] that the L°° estimates also hold for x < 0 
since we can also find the invariant regions in the lower half-plane of u — v, which are 
symmetric with the invariant regions in the upper half-plane around the u axis. 


4. H ^ Compactness 


In this section we shall prove the compactness of the approximate viscous solutions. 
For the strictly convex entropy 

+ 1 


T] = 


and entropy flux 


q = 


L 

-AP + M 


from the parabolic equations ()2.8p . one has 
Vy Qx — ^Xjxx T n(.b, M') 2£ 


L2 


m2 + 1 


L3 ^ I I 


iy 1 (M ~ 

= ^Vxx + n(.b, M) — 2e I -|- y I -yLx ~ Air 

L'* L \ L 



(4.1) 
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where 


From 


we have 


and thus 


n(L,M) 


+ 1 /~ 0 ~ ~ ~0 ~\ 
=^^(rl,L-2rl,M + Ti,N) 

+ ^ - 2fI 2 M + fliiv) . 



V 


M = --, 

V 


L 


X 



Mx 



V 


Vy + Qx 


^Vxx + n(F, M) — 2e 



(4.2) 


Due to the uniform estimates for u and v in Subsections o and 13.21 we have 


0 < bi{5) <L< b2{d), \M\ < b3{6), (4.3) 

uniformly in £ in D, where bi{6), b2{5) and b^{5) are positive constants depending on 5 > 0. 
Therefore n(L, M) is also uniformly bounded in e. Let 


= {(x, 2/) : X G M, - 2/0 <y <0}, 


where 220 > 0 is arbitrary. Choose the test function (j) G C^{Q.) satisfying (()|x = 1, 
0 < <(> < 1, where iL is a compact set and K d S = supp (p. From (j4.2|) . we have 


r r 2£(^+^) <pdydx 

J-00 J-yo \V V J 

< / Uvxx -r]y-qx+ n(L, M) j (pdydx 
J-00 J-yo ^ 

= \ / (£7](j)xx + 'n4>y + Q4>x + 'n{L,M)(j)]dydx 

J —00 J —yo ^ ' 

< M{(p) 


(4.4) 


for some positive constant M{<p) uniform in e G (0,1). Since v is uniformly bounded from 
below with positive lower bound, eu^ are bounded in L)^^(D). Since u is uniformly 
bounded and v has uniform positive lower bound, one has 


eLl = £^< Cevl 


eM^ = e 


UVx 


u. 


< C£{vl + ul), 


for some positive constant C uniform in e, then we see that and sM^ uniformly 
bounded in L)^^(D). Noting that 
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and for arbitrary cj) G 

/ oo rO POO rO 

/ sLxxHydx = / eLx(l)xdydx 

-oo J —vn J —oo J—yo 


' -oo J -yo 

< 

^ C yfs — y 0 QiS s —y 0 , 


eLldydxj ( / {(f)xfdydx 

V^supp (j) J 


(4.5) 


we see that eLxx is compact in From (|2.8I) . we have 

Mx — Ly = f' 22 ^ ~ 2 f^ 2 -^ + fX 1 -/V — eLxx- (4-6) 

Since f' 22 '^ ~ 2 TI 2 M + rf^iV is uniformly bounded in Q, it is uniformly bounded in 
and compact in with some 1 < p < 2 by the imbedding theorem and the 

Schauder theorem. Therefore Mx — Ly is compact in Moreover, we see that 

Mx — Ly is uniformly bounded in (Ll) since M and L are uniformly bounded. 

Finally, by Lemma [4.II below, we conclude that Mx — Ly is compact in Similarly, 

Nx — My is also compact in Since 7 is C^, we see that Mx — Ly and Nx — My 

are also compact in 

We record the following useful lemma (see mm) here: 

Lemma 4.1. Let Ll G M"" be a open set, then (compact set of ) f) (bounded set 

ofWi~^'^{Ll) ) C (compact set of where q and r are constants, 1 < q <2 < r. 


5. Main Theorems and Proofs 


In this section, we shall state our main results and also give the proof. 

In the previous sections, we have established the L°° uniform estimate and com¬ 
pactness of the viscous approximate solutions to (| 2 . 8 p for some special metrics of the 
form 

ds^ = Edx^ + Gdy"^. 

The corresponding Gauss curvature K{x,y) = K has the following form (see [27]): 


1 + ExGx Gx + EyGy 

4.EG y E Gy 


‘2{Eyy -|- Gxx) 


(5.1) 


To prove the existence of isometric immersion, first let us recall the following compensated 
compactness framework in Theorem 4.1 of [ 8 |: 


Lemma 5.1. Let a sequence of functions {L^^, M^, N^){x,y), defined on an open subset 
n C satisfy the following framework: 

(W.l) {L^,M‘^,N‘^){x,y) is uniformly bounded almost everywhere in fl C with respect 
to e; 

(W.2) Mf — Ly and Nf — My are compact in 77^“^ (0); 
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(W.3) There exist o^(l); 3 = 1, 2, 3, with Oj(l) 0 in the sense of distributions as e ^ 0 
such that 

Mf -Ll = Tl^L^ - 2Tl^M^ + + o!(l), 

Nf - M! = + 2T\^M^ - + o|(l), 

and 

= K+ ol{l). (5.3) 

Then there exists a subsequence (still labeled) (L^, , N^) converging weak-star in L°° to 

{L, M, N){x, y) as e —)• 0 such that 

(1) {L,M,N) is also bounded in id G 

(2) the Gauss equation 15.31] is weakly continuous with respect to the subsequence 

(L^, , N^) converging weak-star in L°° to {L, M, N){x,y) as e ^ 0; 

(3) The Codazzi equations 115.^} as e —)• 0 hold for (L, M, N) in the sense of distribu¬ 
tion. 

We now present the results on the existence of isometric immersion of surfaces with the 
two types of metrics studied in Section 3. 

Definition 5.1. A Riemannian metric on a two-dimensional manifold is called a catenoid 
metric if it is of the form 

ds^ = E{y)dx^ + cE{y)dy^, 

with 

c > 0, E{y) > 0, E{y)' < 0 for y < 0, 

and the corresponding Gauss curvature is of the form 

K{y) = -koE{y)-P" 
with constants fco > 0 and ft > \f2. 

From the Definition 15.II and the formula in ()5.ip . we see that E{y) satisfies the following 
ordinary differential equation: 

[E{y)'f - EE{y)" = -2koE{yf-P^. (5.4) 

We can solve it through the following process utilizing the method of [l 2 ]. Set 

E{y) = e^^y\ 

then 

E{yy = e^^y^w{y)y E{y)" = ((u;(y)')' + w{y)") , 

and (j5.4p becomes 

w{y)" = 2koe-^"^^y\ (5.5) 

Let f{w{y)) = w{yy. After differentiating respect to y, we get 

f{wyw{yy = w{yy\ 

i.e., 

f{w)f{wy = 2kQe~^^^, 
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therefore, noting that E' implies w' < 0, 

Then one has 

^ [ dw 

y = C2- / , ^=, 

j ^Ci - 

Denote the right side of the above equation by h{w), then w{y) = h~^{y) and E{y) = 
where h~^{y) is the inverse function of h{w), Ci and C 2 depend on the the value 
of w(0) and w(0)'. Then the catenoid metric is of the form: 

ds^ = + ce^ (5.6) 

We note that the special catenoid type metric given in Example 13.11 

E(y) = (ccosh(^))^, G{y) = ^^^^-^E{y), 

with K{y) = — l)£'(y)“^^ is a catenoid metric in the sense of Definition 15.II 

Definition 5.2. A Riemannian metric on a two-dimensional manifold is called a helicoid 
metric if it is of the form 

ds'^ = E{y)dx^ + dy^, 

with 

E{y) > 0, E{yy < 0 for 1/ < 0, 
and the corresponding Gauss curvature is of the form 

K{y) = -koE{yf 

with constants A:o > 0 and a < —2. 


From Definition 15.21 and the formula (j5.ip . E{y) satisfies the following ordinary differ¬ 
ential equation: 

{EiyYf - 2EE{y)" = -AkoE{yf+\ (5.7) 


Set 


E{y) = w{yf, 


then (15.71) becomes 


w{y)" = kow{yf<^+\ 

Letting g{w{y)) = w{y)', and differentiating respect to y, one has 

w{yf = u{w)'w{y)\ 


i.e., 

Thus 


g{w)'g{w) = kow 


2a+l 


dw 


= Ci + 


u -|- 1 


-w 


2a+2 
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Then 


y = C2- 


dw 




Denote the right side of the above equation by h{w)^ then w{y) 
Therefore, the helicoid metric is 

ds^ = {h~^ {y))‘^ dx^ + 

where h~^{y) is the inverse function of 


h{w) = C 2 


dw 


h i(y), and E{y) = 


Similar to the catenoid metric, Ci and C 2 depend on the value of r(;(0) and 'w{0y. As an 
example, the helicoid surface with 

E{y) = c^ + y^, K{y) = - ^ 

+ y-^Y 

is a two-dimensional Riemannian manifold with the helicoid metric in the sense of Dehni- 
tion 15.21 

Now we prove the existence of isometric immersion of surfaces with the above two 
types of metrics into by using Lemma l5.ll 


Theorem 5.1. For any given yo > 0, let the initial data 

(u, v)\y=-y^ = (ho(x), i;o(x)) := (n(x, -yo),v{x, -yo)) (5.8) 

satisfy the following conditions: 

uq + vq and uq — uq are bounded, 


and 

inf (ho + ho) > 0, sup (ho — ho) < 0. 

XeR 3;gR 

Then, for the catenoid metric in the sense of Definition 15.ll the Gauss-Codazzi system 
112.1\) has a weak solution in Q = {(x,y) : x G M, —yo < y < 0} with the initial data 15.dl) . 


Remark 5.1. For the initial data (no(x), uo(x)) satisfying the conditions in Theorem 15.11 
there exists a constant 6 > 0, such that 

(ho(x), ho(x)) G ACBD := {{u, v) : 6<w = u + v< 2uo + 6, 

— {2uq + 5)<z = u — v< —h}, 

for the catenoid type metrics, where uq is dehned in ()3.12p . That is, (ho(x), ho(x)) lies in 
the invariant region ACBD sketched in Figure 3 or Figure 5, i.e., ho(x) is bounded, and 
ho(x) has positive lower bound and upper bound. 


Proof First for the initial data (|5.8p . the corresponding initial data for L,M,N is 


Lo{x) = - 


1 


Ar„(.) = 


No{x) = 


UoixY - Vo{x) 


Uo(x)’ Voix) ’ vo{x) 

We use system (|2.19p to obtain the approximate viscous solutions and their L°° estimate, 
and then use (j2.8p to obtain the compactness. 
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Step 1. We mollify the initial data (|5.8p as 

•uola:) = uo{x) * f, Vo{x) = voix) * f, 

where is the standard mollifier and * is for the convolution. Prom the above Remark 
EH there exists a (5 > 0, such that, 

1110(^)1 — 0 < (5 < Vq{x) < M{6), 

where M{6) is positive constant depending on 5, and 

(no(2^)) f’o(®)) ^ {uo{x),vo{x)) as e ^ 0, a.e. 


Step 2. The local existence of (12.8p can be obtained by the standard theory, hence 
we have the local existence of (|2.19l) . From Section 3, we have the L°° estimate for the 
approximate solution of (|2.19l) , then we can obtain the global existence of the approximate 
solution as follows. We observe that the first equation of (I2.19P is of the divergence form, 
thus the estimate of can be handled in the standard way. However, the second equation 
of (I2.19P is not in divergence form. By differentiation of the second equation with respect 
to X, we get 

iyx)y T {uVx VUx)x — T 1 ; 

where we omit e and still use u{x,y),v{x,y) for the approximate solution of (|2.19l) . Let 
G{x,y) is the heat kernel of hy = ehxx, then we can solve the above equation as the 
following: 

/ OO 

G{x - ^,y)voxiOd^ 

-OO 

+ j J (^9{u,v)x - {uvx-vux)x + ^ G{x - ^,y - r])dCdr] 


G{x - ^,y)v^ 0 xiOd^ 


+ 


f I 

'J —yo 'J —OO 


-9{u, v) + {uVx - VUx) + 


2evt, 


G{x -i,y - r])xdCdr]. 


Therefore 


< 


/ OD 

G{x - ^,y)\vl,x{C)\d^ 

-OO 

fO poo / 2fqp \ 

+ J J i^\ 9 {u,v)\ + \{uvx-vux)\^ - ^ j G{x - ^,y -r])x\d^dr] 

/ O /•OO 

/ \Gix - C,y - r])x\d^dr], 

-yo J -OO 


where || • 11(^0 stands for the norm of continuous functions. The above process yields 
the G^ estimate of Vx- Similarly, we can estimate G^ (for any integer A: > 1) norms of 
u{x,y),v{x,y), which can be bounded by the norms of Uq and Vq. Thus, we obtain the 
global existence of smooth solutions to the system (j2.19p . 

Step 3. In Section 3 we have proved boundness of L,M,N, and then L, M, N is 
in L°° for 7 £ By the reverse process of Section 2, we can reformulate the equations 
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(12.191) of u and v as the equations (|2.8I) . Therefore, as in Section 4 we also obtain the 
compactness. So we have proved that our approximate solutions satisfy (W.l) and (W.2) 
in the framework of Lemma l5.II Furthermore, from Section 2, we have 

(M, - Ly) - - 2 tl^M + 

{N, - My) - + 2 r\^M - fJiiV) = 

As in Section 4, from (14. Sp . we can get that, as e —>• 0, 

0 , 

in the sense of distribution, and then 

Mx - Ly = tl^L - 2V\2M + + oi(l) 

holds in the sense of distribution. Similarly, 

Nx - My = -f^ 2 ^ + 2 f laM - r\iN + 02(1) 

also holds in the sense of distribution. Here oi(l), 02 (l) —>■ 0 as e ^ 0. We note that the 
Gauss equation holds exactly for the viscous approximate solutions. Therefore (W.3) is 
satisfied. Consequently, we complete the proof of the theorem and obtain the isometric 
immersion of the surface with the catenoid metric in using Lemma l5.II □ 


Since we also obtained the L°° estimate and the compactness for the helicoid 
metric in the previous sections, we can have the isometric immersion of surfaces with the 
helicoid metric in just as the case for the catenoid metric. 

Theorem 5.2. For any given yo > 0, let the initial data 

{u,v)\y=-yo = {uo{x),vo{x)) := {u{x,-yo),v{x,-yo)) 
satisfy the following conditions: 

uq + vq and uq — uq are bounded, 

and 

inf (no + Vq) > 0 , sup (no — vo) < 0 . 

3:GR a:GR 

Then for the helicoid metric in the sense of Definition \5. 21 the Gauss-Codazzi system 112.1\) 
has a weak solution in kl = {{x,y) : x G R, —yo < 2 / < 0} with the initial data L5.8]) . 

Remark 5.2. Although the catenoid metric for (3 > y/2 has also been studied by Chen- 
Slemrod-Wang in [ 8 ], their estimate is different from ours, especially for L. In addition, 
we also prove the isometric immersion of catenoid for [3 = y/2. 

Remark 5.3. By Remark 13.11 if we replace the second condition on the initial data in 
Theorem 15.11 and Theorem 15.21 by the following: 

sup (no + flo) < 0, inf (no - flo) > 0, 
xGR *6® 

then both theorems still hold. 

Remark 5.4. By even symmetry from the weak solution in H = {(x,y) : x G R, —yo < 
y < 0} obtained in Theorem 15.11 and Theorem 15.21 we can obtain the weak solution in 
^2' = {{x,y) : X G R, 0 < y < yo}. 
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